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ABSTRACT 

We report on a ^12 hr XMM-Newton observation of the supergiant High-Mass X-ray Binary IGR J 16207- 
5 129. This is only the second soft X-ray (0.4-15 keV, in this case) study of the source since it was discovered by 
the INTEGRAL satellite. The average energy spectrum is very similar to those of neutron star HMXBs, being 
dominated by a highly absorbed power-law component with a photon index of F = LlSiliogj. spectrum 
also exhibits a soft excess below ^2 keV and an iron Ka emission line at 6.39 ±0.03 keV. For the primary 
power-law component, the column density is (1.19;!;Q g5) ^ ^^^^ cm"^, indicating local absorption, likely from 
the stellar wind, and placing IGR J16207-5129 in the category of obscured IGR HMXBs. The source exhibits 
a very high level of variability with an rms noise level of 64%±21% in the lO"** to 0.05 Hz frequency range. 
Although the energy spectrum suggests that the system may harbor a neutron star, no pulsations are detected 
with a 90% confidence upper limit of ^2% in a frequency range from ^lO"** to 88 Hz. We discuss similarities 
between IGR J16207-5129 and other apparendy non-pulsating HMXBs, including other IGR HMXBs as well 
as 4U 2206H-54 and 4U 1700-377. 

Subject headings: stars: neutron — black hole physics — X-rays: stars — stars: supergiants — stars: individual 
(IGR J 16207-5 129) 



1. INTRODUCTION 

The hard X-ray imaging by the INTEGRAL satellite 
(IWinkler et al.ll20()3h has been uncovering a large number of 
hard X-ray sources. Since INTEGRAL'S launch in 2002 Oc- 
tober, 550 sources have been detected by the IBIS instru- 
ment in the ^20-50 keV band (based on version 29 of the 
"General Reference Catalog"). Included in these sources are 
236 "IGR" sources that were unknown or at least not well- 
studied prior to INTEGRAL. An important result from the IN- 
TEGRAL mission has been the discovery of a relatively large 
number of High-Mass X-ray Binaries (HMXBs). There are 
37 IGR sources that have been classified as HMXBs (and it 
should be noted that about 1/3 of the IGR sources are still 
unclassified). The IGR HMXBs are interesting both for the 
large number of new systems as well as the specific prop- 
erties of these systems. These include a new chss o f "Su- 
pergiant Fast X-ray Transients" (iNegueruela et alJl2006,) that 
are HMXBs that can exhibit hard X-ray flares that only last 
for a few hours while the X-ray flux changes by orders of 
magnitude (in 't Zand 2005; Sguera et al. 2006). Many of the 
IGR HMXBs are also extreme in having a high level of ob- 
scuratio n (A^h^'^IO^'^t-IO^'^ crn~^) due to mater ial local to the 
source dWalter et all 120061: IChatv et all [20081) . For both the 
SFXTs and the obscured HMXBs, it is thought that a strong 
stellar wind is at least partially res ponsi ble for their extreme 
X-ray prop erties (Filliatr e & Chatvll2004t (Walter et alj|2006l: 
IWalter & Zurita Heras,.2007l) ! 
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The 37 IGR HMXBs presumably contain either neutron 
stars or black holes, but the nature of the compact object is 
only clear for the 12 systems for which X-ray pulsations from 
their neutron stars have been detected. Eleven of these sys- 
tems have pulse periods ranging from 5 to 1300 s, a nd the 12th 
system has an unusually long period of ^5900 s dPatel et al] 
2004). Another X-ray property that can be taken as evidence 
for the presence of a neutron star is a very hard X-ray spec- 
trum. Typically, neutron star HMXBs have X-ray spectra that 
can be modeled with a power-law with a photo n index of F'^ 1 
that is exponentially cutoff near 10-20 keV ('NagasdliasJ; 
Lutovinov et al. 2005). Although there are only a few known 
black hole HMXBs (e.g., Cygnus X-1, LMC X-1, LMC X- 
3, M33 X-7), their X-ray spectral properties are similar to 
the general class of black hole X-ray binaries with power- 
law spectra with F^ 1.4-2.1 in their hardest spectral state 
(McClintock & Remillard 2006). If black hole spectra show 
a cutoff, it is usuaUy close to 100 keV (iGrove et alJ ll998). 
significantly higher than seen for neutron star HMXBs. How- 
ever, it should be noted that evidence from the shape of the 
spectral continuum alone is usually taken only as an indi- 
cation of the nature of the compact object. For example, 
HMXBs 4U 1700-37 and 4U 2206H-54 both have continuum 
X-ray spectra similar to neutron star HMXBs, but they are 
considered to be only probable neutron star systems because 
pulsations have not been detected. 

In this study, we focus on X-ray observations of 
IGR J 16207-5 129, which was discovered in the Norma re- 
gion of the Galaxy relat i vely early-on in the I NTEGRAL mis- 
sion (IWalter et alj|2004t iTomsick et alj|2004 . Alfliough it is 
a relatively faint hard X-ray source at 3.3 ±0.1 millicrab in 
the 20^0 keV band dBird et alj 120071) . it has been consis- 
tently detected by INTEGRAL as well as in X-ray follow-up 
observations by Chandra and XMM-Newton (this work), in- 
dicating that it is a persistent source. The Chandra observa- 
tion provided a sub-arcsecond position that allowed for the 
identification of an optical counterpart with R = 15.38 ±0.03 
and an IR counterpart with Ks = 9.13 ± 0.02 dTomsick et alj 
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12006). Based on the op tical/IR Spectral Energy Distribu- 
tion. iTomsick et all (120061) found that the system has a mas- 
sive O- or B-type optical companion. Optical and IR obser- 
vations confirme d this and indicat e a su pergiant nature for 
the companion (Masetti et al. 2006; Neaueruela & Schurch| 
l2007HRahouiet al . 2008). and Rahoui et al. (2008) estimated 
a source distance of ~4.1 kpc. With further IR spectroscopy, 
the spectral type of the companion was narrowed down to 
Bl la and a source distance of 6.1;!;^ '^ kpc was estimated 
dNespoli. Fabregat & Menn ickerit 2008). 

Our current knowledge about the soft X-ray properties of 
IGR J 16207-5 129 comes from the Chandra observation that 
was made in 2005. This showed that the source has a hard 
X-ray spectrum with a power-law photon index of F = 0.5^° 5, 
and the source also exhibited significant variability over the 
5 ks Chandra observation (Tomsick et al. 2006). Based on 
these properties, we selected this target for follow-up XMM- 
Newton observations to obtain an improved X-ray spectrum 
and to search for pulsations that would provide information 
about the nature of the compact object. Here, we present the 
results of the XMM-Newton observation. 

2. XMM-NEWTON OBSERVATIONS AND LIGHT CURVE 

We observed IGR J 16207-5 129 with XMM-Newton dur- 
ing satellite revolution 1329. The observation (ObslD 
0402920201) started on 2007 March 13, 8.27 hr UT and lasted 
for44ks. The EPlC/pn instrument ( Striider et al. 2001)accu- 
mulated ^0.4-15 keV photons in "small window" mode, giv- 
ing a 4.4-by-4.4 arcminute^ field-of-view (FOV) and a time 
resolution of 5.6718 ms. The mode used for the 2 EPIC/MOS 
units (Turner et al. 2001) is also called "small window" mode, 
and for MOS, its features are a 1.8-by-1.8 arcminute^ FOV 
and a time resolution of 0.3 s. In addition to the XMM-Newton 
data, we downloaded the "current calibration files" indicated 
as necessary for this observation by the on-line software tool 
cif build. For further analysis of the data, we used the 
XMM-Newton Science Analysis Software (SAS-8.0.0) as well 
as the XSPEC package for spectral analysis and IDL for tim- 
ing analysis. 

We began by using SAS to produce pn and MOS images 
and found a strong X-ray source consistent w ith the Chandra 
position of IGR J16207-5129 ( ITomsick et aLll2006h . This is 
the only source seen in the pn image, but for MOS, five of the 
outer CCD detectors are active, and a few other faint sources 
are detected. To produce a pn light curve, we used SAS to 
read the event list produced by the standard data pipeline. 
We extracted source counts from a circular region with a ra- 
dius of ^35", which includes nearly all of the counts from 
the source. For subtracting the background contribution, we 
used the counts from a rectangular region with an area 3.5 
times larger than the source region located so that no point 
in the background region comes within 2' of the source. In 
producing the light curves, we applied the standard filtering 
("FLAG=0" and "PATTERN <4") as well as restricting the 
energy range to 0.4-15 keV. 

The pn light curves with 50 s time resolution are shown in 
Figure [T] While the average count rate after deadtime correc- 
tion and background subtraction is 1.64 c/s, there is a very 
high level of variability with count rates in the 50 s time bins 
ranging from 0.03±0.15c/sto7.6±0.5c/s. The background 
light curve after deadtime correction and scaling to the size of 
the source region is shown in Figure[T}5. The background rate 
also shows a high level of variability during the observation. 
The average background rate for the entire 44 ks observation 
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Fig. 1 .— (a) XMM-Newton pn 0.4^15 keV light curve for IGR JI6207- 
5129. The time resolution is 50 s, and we have subtracted the background 
contribution, (b) The light curve from the background region after scaling 
to the size of the source extraction region. The solid hne marks the low- 
background time segment used for spectral analysis. 

is 0.40 c/s, but it is much higher at the end of the observation. 
For the first 39 ks, the average background rate is 0.18 c/s, 
while it is 1.9 c/s for the last 5 ks. We verified that the higher 
count rate is due to proton flares by producing a pn light curve 
for the full pn FOV in the >10 keV energy band. The average 
full-FOV count rate in this energy band is 0.33 c/s for the first 
39 ks while it is 4.9 c/s for the last 5 ks. Thus, for the spectral 
analysis described below, we only included the portion of the 
exposure indicated in Figure [T] 

We also produced MOS light curves using a circular source 
region with a radius of ~35", and they show the same vari- 
ability and flaring as the pn light curves. We used the MOSl 
and MOS2 pipeline event lists, and applied the standard fil- 
tering ("FLAG=0" and "PATTERN<12"). In small window 
mode, the active area of the central CCD is too small to use 
the data from this CCD for background subtraction; thus, we 
used a source-free rectangular region from one of the outer 
CCDs as our background region. We also made MOS light 
curves using data from the full FOVs in the >10 keV energy 
band. They show that a large increase in the background level 
occurred simultaneously with the increase seen in the pn, and 
for the analysis described below, we used the MOS data from 
the first 39 ks of the observation (i.e., the low -background 
time indicated in Figure [TJ)). After background subtraction, 
the 0.4-12 keV MOSl and MOS2 count rates in the source 
region during the low -background time are 0.50 and 0.51 c/s, 
respectively. 

The third X-ray instrument on XMM-Newton is the Reflec- 
tion Grating Spectrometer (RGS). We inspected the RGS dis- 
persion images that are used to extract spectra, but there is no 
evidence that IGR J 16207-5 129 is detected. To determine if 
an RGS upper limit is constraining, we used the spectra ob- 
tained using the EPIC instruments (see §3.2 below) along with 
an RGS response matrix. We find that, for the entire 44 ks ob- 
servation, we would expect the RGS to collect ^22 counts and 
^13 counts in the first and second grating orders, respectively, 
which, given the instrumental background, is consistent with 
the non-detection. 



Obscured HMXB IGR J16207-5129 



3 



3. ANALYSIS AND RESULTS 
3.L Timing 

We used the XMM-Newton instrument with the highest ef- 
fective area, the pn, for timing analysis and started with the 
standard pipeHne event Hst. We used the SAS tool barycen 
to correct the timestamps to the Earth's barycenter for each 
event. As IGR J 16207-5 129 is a HMXB that may contain a 
pulsar, a primary goal of this analysis is to search for peri- 
odic signals. We used the IDL software package to read in 
the event list and make a 0.4-15 keV light curve at the high- 
est possible time resolution, Af = 5.6718 ms. We note that it 
is important to use this exact value for Af to avoid producing 
artifacts in the power spectrum. Once the high time resolution 
light curve was produced, we used IDL's Fast Fourier Trans- 
form (FFT ) algorithm to pro duce a Leahy-normalized power 
spectrum ( Leahy et al.lll983h . The power spectrum extends 
from 2.3 x 10"^ Hz (based on the 44 ks duration of the obser- 
vation) to 88 Hz (the Nyquist frequency). 

Part of this power spectrum is shown in Figure |2] The 
power is distributed as a probability distribution with 2 
degrees of freedom (dof), and we used this fact along with the 
total number of trials, which is equal to the number of fre- 
quency bins (Niriais = 3,839,999 in this case) to determine a 
detection threshold. Here, the 90% confidence detection limit 
is at a Leahy Power of 34.9, and this is shown in Figure |2] 
This limit is not exceeded at frequencies above 0.005 Hz, but 
there are a large number of frequency bins below 0.005 Hz 
that exceed the limit. We suspect that this may be related to 
low-frequency continuum noise; thus, we treat the >0.005 Hz 
case first. 

The highest power that we measure in the 0.005-88 Hz 
range is 33.4, which is just below the 90% confidence detec- 
tion threshold. We do not consider this as even a marginal 
detection, but we use this value (Pmax) to calculate an up- 
pe r limit on the streng th of a periodic signal. As described 
in Ivan der the 90% confidence upper limit is 

given by Pul = Pmax-Pexceed, whcrc Pexceed IS the powcr level 
that is exceeded in 90% of the frequency bins. In our case, 
Pexceed = 0.2 SO that Pm = 33.2, which, after converting to 
fractional rms units using the average source and background 
count rates, corresponds to an upper limit on the rms noise 
level for a periodic signal of <2.3%. We also produced a 
power spectrum with 4.6 x lO"*" Hz frequency bins, but we 
still did not find any bins with power in the 0.005-88 Hz range 
with powers exceeding the 90% confidence upper limit. 

To characterize the low-frequency noise, we produced a 
0.4-15 keV pn light curve with 10 s time resolution, and made 
a new power spectrum with a Nyquist frequency of 0.05 Hz. 
The new power spectrum consists of the average of power 
spectra from four ^ 1 1 ks segments of the light curve, giv- 
ing a minimum frequency of 9.2 x 10"^ Hz. In addition, we 
converted the power spectrum to rms normalization and re- 
binned the power spectrum so that the power in each bin fol- 
lows a Gaussian distribution, allowing for us to fit a model 
to the power spectrum with ^^-minimization. The resulting 
power spectrum is shown in Figure |3] and the figure also 
shows that the power spectrum is well-described by a power- 
law model (P = A{i'/l Hz)~"). We obtain a power-law in- 
dex of a = 1.76 ±0.05, and the integrated rms noise level is 
64%±21% (0.000092-0.05 Hz). For the fit, the = 43.6 for 
39 dof, indicating that the power-law provides an acceptable 
description of the power spectrum. 

To search for a periodic signal in the presence of red 
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Fig. 2. — The lowest frequency portion of the Leahy-normahzed power 
spectrum for IGR J16207-5129 with a frequency binsize of 2.3 X 10"^ Hz. 
The horizontal dashed line marks the 90% confidence detection limit (after 
accounting for trials). The full power spectrum extends to the Nyquist fre- 
quency of 88 Hz, based on the pn time resolution of 5.6718 ms, and only 
exceeds the 90% confidence detection Umit in the region below 0.05 Hz. 
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Fig. 3. — Red noise rms-normalized power spectrum for IGR J16207- 
5129. The power-spectrum is fitted with a power-law (solid Hne). 

(i.e., power-law) noise requires an extra step in the analy- 
sis. Specifically, after producing another Leahy-normalized 
power spectrum, this time covering the 0.000023-0.05 Hz 
frequency range, multiplying this power spectrum by 2 and 
dividing by the (re-normalized) power-law model leads to a 
power spectrum where the power i n each frequency b in fol- 
lows a x^ distribution with 2 dof dvan der Klislll989l) . This 
is illustrated by showing these power spectra in Figure]?] and 
we can now search the powers in the bottom panel for periodic 
signals. After accounting for the number of trials (i.e., 2,176 
frequency bins), the 90% confidence detection limit is shown. 
Although there are no signals that reach this detection limit, 
the maximum signal at 0.0089 Hz (112.56 s) is only slightly 
below. We do not consider this to be even a marginal detec- 
tion, but we use the power in this bin, P„ax = 19.2 to determine 
the upper limit. Using the same procedure as described above 
results in an upper limit on the rms noise level for a periodic 
signal of <1.7%. 

One caveat to the non-detection of a periodic signal is that 
it is possible for the power from a periodic signal to be spread 
out in frequency by orbital motion. This can happen if the du- 
ration of the observation is a substantial fraction of the orbital 
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Fig. 4. — (a) Leahy-normalized power spectrum for IGR J16207-5129 with 
a frequency binsize of 2.3 X 10"^ Hz. The solid line is the power-law model 
from the fit shown in Figure|3] (b) Two times the data-to-model ratio using the 
data and model shown in panel a. The horizontal dashed line shows the 90% 
confidence detection limit (after accounting for the number of trials). The 
strongest signal is at 0.0089 Hz (1 12.56 s). but it is not statistically significant. 



Fig. 5. — XMM-Newton pn and MOS energy spectrum for IGR J16207- 
5129 fitted with an absorbed power-law. Residuals (lower panel) illustrate 
the presence of a soft excess. The pn points are marked with points while 
the MOSl and MOS2 points are not (upper and lower panels), and the line 
used for the pn model is dashed while the MOS model lines are solid (upper 
panel). 



period. For IGR J 16207-5 129, we do not know the orbital 
period; however, we do know that it is an HMXB and is ex- 
pected to have an orbital period of between several days and 
several weeks. Thus, our ^^12 hr XMM-Newton observation 
should cover only a small fraction of the orbit. 

3.2. Energy Spectrum 
3.2.1. Time-Averaged Spectrum 

We extracted pn, MOSl, and MOS2 energy spectra using 
the source and background regions and filtering criteria de- 
scribed above, and we produced response matrices using the 
SAS tools rmf gen and arf gen. We obtained a total expo- 
sure time of 27 ks for the pn and 37 ks for each MOS unit. 
The exposure time is lower for the pn due to a higher level 
of deadtime. We rebinned the energy spectra by requiring at 
least 50 counts for each energy bin, leaving a pn spectrum 
with 722 bins, a MOSl spectrum with 290 bins, and a MOS2 
spectrum with 298 bins. Although we use the 0.4-15 keV en- 
ergy band for pn light curves, we restrict the spectral analysis 
to 0.4-12 keV for the pn and MOS units as the caUbration^ 
has focused on this energy range. 

We used the XSPEC software package to jointly fit the 
3 spectra with an absorbed power-law model. To account 
for absorption, we used the photoelectric abso rption cross 
sections from B alucinska- Church & McCammon (1992) and 
elemental abundances from Wilms, Allen & McCray (2000,). 
which correspond to the estimated abundances for the inter- 
stellar medium. We left the relative normalization between 
the instruments as free parameters, and a x^-minimization 
fit indicates that the MOS normalizations are slightly higher 
(3%±1%) than the pn normalization. The absorbed power- 
law model, with the best fit values of A^h = 9 x 10^^ cm"^ for 
the column density and F = 0.97 for the power-law photon 
index, appears to provide a good description of the spectral 
continuum above 2 keV; however, the fit is not statistically 
acceptable overall with X' lv= 1778/1305. The counts spec- 
trum and residuals for the absorbed power-law fit are shown 

* A document on the pn and EPIC calibration 
dated 2008 April by M. Guainazzi can be found at 
|http://xmm2.esac.esa.uit/docs/documents/CAL-TN-0018.pdf| 
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Fig. 6. — XMM-Newton pn and MOS energy spectrum in the iron line 
region for IGR J 16207-5 129 fitted with an absorbed power-law. Residuals 
(lower panel) illustrate the presence of an iron line. The pn points are marked 
with points while the MOS 1 and MOS2 points are not (upper and lower pan- 
els), and the line used for the pn model is dashed while the MOS model lines 
are solid (upper panel). 

in Figure|5] Large residuals are seen below 2 keV, suggesting 
the presence of a soft excess. We also produced pn and MOS 
spectra with a higher level of binning to study the iron Ka 
region of the spectrum. The spectrum and residuals shown in 
Figure |6] indicates the presence of an iron emission line near 
6.4 keV. 

We refit the spectrum (going back to the lower level of 
binning) after adding a Gaussian to model the iron line and 
a second spectral component to account for the soft excess. 
We tried different models for the second spectral compo- 
nent, including a Bremsstrahlung model, a black-body, and 
a power-law. In each case, the second component was ab- 
sorbed, but we fixed A^ h for this component to th e Galactic 
value, 1.7 x 10^^ cm'^ (iDickey & Lockmanl[T990h . We also 
used this level of absorption for the iron line. Adding the iron 
line and using a Bremsstrahlung model for the second com- 
ponent yields a much improved fit (over the power-law alone) 
of x^/^ = 1402/1300. The best fit Bremsstrahlung temper- 
ature is 189 keV, which is approaching the upper end of the 
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Fig. 7. — XMM-Newton pn and MOS unfolded energy spectrum for 
IGR J 16207-5 129 with a model consisting of an absorbed power-law (dashed 
line), a less absorbed power-law to account for the soft excess (dotted line), 
and an iron K« emission line (solid line). 

range allowed by the XSPEC model bremss. We derive a 
90% confidence lower limit on the temperature of 50 keV, in- 
dicating that if the soft excess is due to Bremsstrahlung emis- 
sion, it is rather hot. While it is tempting to take the high 
Bremsstrahlung temperature as an indication that the soft ex- 
cess is non-thermal, using a black-body model for the sec- 
ond component yields a fit of comparable quality, 'X^ jv = 
1403/1300, and the black-body temperature is ^^0.6 keV. Fi- 
nally, if a power-law is used for the second component, the 
quality of the fit is identical to the Bremsstrahlung fit, x^ jv = 
1402/1300, and the power-law photon index is T = 0.9!|] j. 
Although the statistical quality of the fits does not allow us to 
determine which model is the best to use for the second, soft 
excess, component, the fact that the power-law has a photon 
index that is consistent with the value of F found for the pri- 
mary power-law component, F = 1.15;!;j]'Q5 (90% confidence 
errors), suggests that it may be possible to interpret the soft 
excess as an unabsorbed portion of the primary component, 
and this possibility is explored further in §3.2.2. Thus, we 
proceed by focusing on the model that uses the power-law for 
the second component. 

The different components of the model with an absorbed 
power-law, a power-law with just interstellar absorption, and 
an iron line are shown in Figure |7] and the model parameters 
are shown in Table [T] The primary power-law component has 
the value of F given above absorbed by a column density of 
A^H = (l-19!o:!]5) X 10^^ cm-2, and the 0.5-10 keV unabsorbed 
flux of the power-law component is 3.7 x 10"'' ergs cm~^ s~', 
which is a factor of ^-^20 higher than the soft excess compo- 
nent. The parameters of the narrow iron Ka line are well- 
constrained. The energy of the line is 6.39 ±0.03 keV, which 
is consistent wi th iron ion ization states between neutral (Fel) 
and ^FeX (INagase et al.lll986l) . The 90% confidence upper 
limit on the width of the line is <0.12 keV. The line is clearly 
detected, but it is not extremely strong with an equivalent 
width of 42 ± 12 eV after correcting the line and power-law 
components for absorption. 

3.2.2. Spectral Properties vs. Intensity 

A caveat on the above spectral analysis is that the parame- 
ters represent average spectral properties for a highly variable 
source. As an initial look at how the spectrum changes with 
flux level, we extracted deadtime corrected and background 
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Fig. 8.— IGR J16207-5129 light curves in 0.4-2 keV (a), 2-5 keV (b), 
and 5-15 keV (c) energy bands and the hardness vs. time (d). Defining the 
5-15 keV rate as Chard and the 2-5 keV rate at C,o/Y, the hardness is (Chard- 
CsoftViCiiard+Csoft)- The time resolution in each case is 100 s, and we have 
subtracted the background contribution. 

subtracted pn light curves with 100 s time resolution in the 
0.4-2 keV, 2-5 keV, and 5-15 keV energy bands. Figure [8] 
shows these light curves (panels a-c) along with the hardness 
vs. time (panel d). The hardness is calculated using the rates in 
the 2-5 keV and 5-15 keV energy bands (defined as Csofi and 
Ch„rd, respectively) according to {Chard-C.,oft)f(Chard+Qoft)- 
The most striking aspect of Figure [8] is how similar the 2- 
5 keV and 5-15 keV light curves appear, suggesting that the 
variability is relatively independent of energy. Although the 
statistics are poorer for the 0.4-2 keV light curve, many of 
the features in this light curve can be associated with variabil- 
ity in the higher energy light curves. A close inspection of 
the hardness (Figure [8]l) shows that the light curves do have 
some energy dependence. While there are some exceptions, 
for most of the deepest dips, the spectra are softer, with the 
hardness being less than zero for many dips. In contrast, most 
flares are harder, with typical hardness values of ^0.2-0.3. 
While the sharp flare that occurs at a time of ■^16,000 s fol- 
lows the typical behavior, a notable exception is found for the 
flare seen in the first ^3,000 s of the observation, which has 
hardness values between -1 and ^0. 

The same hardness values are plotted vs. the 2-15 keV 
count rate (i.e., intensity) in Figure |9l with each point cor- 
responding to a 100 s time bin. In addition, we calculated 
weighted averages in seven intensity bins, and these are also 
shown in the figure. At the high intensity end, the average 
hardness is dominated by the softer flare at the beginning of 
the observation. Then, the spectrum is harder at intermedi- 
ate intensities before softening again below 1 c/s. We cal- 
culated the expected hardness change if the variability is due 
only to a change in A^h, and the prediction as A^h changes from 
1.7 X 10^^ cm"^ to lO^"* cm"^ is shown in Figure|9l The pre- 
dicted change in hardness is drastically larger than the change 
we observe, indicating that changing A^h cannot be the sole 
cause of the variability. 
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Fig. 9. — Hardness-Intensity diagram for IGR J16207-5129. We used 
the 100 s time resolution light curves for the counts in the 5-15 keV band 
(Chard) and the 2-5 keV band (Csoft)- The hardness is defined as (Chard- 
Csoft)/(Ci,ard+Csaft)- In addition to the individual 100 s points, which are 
grey in color, the averages in seven intensity bins are shown (the black his- 
togram). The dashed line shows the hardness-intensity relationship expected 
if the change in count rate was due only to a change in column density. 

To investigate which spectral parameters do change with in- 
tensity, we used the 0.4-15 keV pn light curve to divide the pn 
and MOS data into spectra at four intensity levels: 4-8 pn c/s, 
2-4 pn c/s, 1-2 pn c/s, and 0-1 pn c/s. We used all of the data 
from the first 39 ks of the observation and obtained pn expo- 
sure times of 2,485 s, 6,177 s, 8,520 s, and 9,940 s and MOS 
exposure times of 3,412 s, 8,482 s, 1 1,700 s, and 13,650 s for 
the four intensity levels, respectively. In fitting these spectra, 
we simplified the two power-law model by requiring the same 
r for both power-laws, which is equivalent to using a model 
where a fraction, /, of the power-law component is absorbed 
both by interstellar material and by the extra material local to 
the source, while a fraction, l-f is absorbed only by interstel- 
lar material. Table |2] shows the parameters obtained for the 
total spectrum and also for each of the four intensity levels, 
and Figure [To] shows the unfolded spectra for the four levels. 
While we performed all of the fits with pn and MOS spectra, 
Figure[TO]shows only the pn spectra for clarity. 

One interesting result that we obtain from the spectral pa- 
rameters shown in Table |2] is that, at A^h = (6.4 ±0.4) x 10^^ 
cm"^, the column density for the highest intensity level is 
significantly less that for the lower levels. This spectrum is 
dominated by the flare that occurs during the first 3,000 s of 
the observation, and this indicates that lower A^h is the rea- 
son that this flare is softer. The parameters also indicate that 
there are two reasons why the spectrum is softer during the 
deepest dips. First, the power-law index indicates that the 
source gradually softens toward the lower levels, going from 
r= 1.06 ±0.05 at 4-8 c/s to 1.30 ±0.12 at 0-1 c/s. Secondly, 
the parameters and the spectra shown in Figure [TO] indicate 
that the flux of the soft excess does not decrease as much as 
the flux of the primary power-law component. Between 2-4 
c/s and 0-1 c/s, the flux of the primary power-law component 
changes by a factor of 5. 1 ± 0.5 while the flux of the soft ex- 
cess changes by a factor of 2.0 ± 0.5. The iron line parameters 
show that the equivalent width of the iron line is consistent 
with being constant while the flux of the line changes with 
the continuum flux. 

4. DISCUSSION 
4.1. Is IGR J16207-5129 an Obscured HMXB? 
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Fig. 10. — Unfolded pn energy spectra after dividing the data into the 
following different intensity levels: (a) 4-8 pn c/s; (b) 2-A pn c/s; (c) 1-2 pn 
c/s; and (d) 0-1 pn c/s. In each panel, the best fit partial covering fraction 
models are plotted for all four intensity levels. 

While the Chandra position and the optical/IR spec- 
troscopy leave essentially no doubt that IGR J 16207-5 129 
is an HMXB, it has not been entirely clear whether the col- 
umn density is high enough to require that part of the ab- 
sorbing material is local to the source or not. The only 
other soft X-ray measurement besides the one reported here 
is the Chandra observation, and Tomsick et al. (2006) re- 
ported a value of A^h = (3. 7!':^) x 10^^ ^m-^ (90% confi- 
dence errors), while many of the obscured HMXBs have col- 
umn densities in excess of 10^^ cm"^. When one considers 
atomic and molecular hydrogen, the total Galactic column 
density is near 2.4 x 10^ ' cm"^ (Di ckey & Lockman 199^ 
iDame. Hartmann & Thadd eus 2001: ITomsick et al.ll2008l) . so 
the Chandra measurement is only marginally higher than the 
Galactic value. 

However, the value that we measure with XMM-Newton 
is substantiafly higher, A^h = (l-19;l:oo5) ^ ^^^^ cm"^, which 
is significantly in excess of the Galactic value. To 
check on whether the difference between the Chandra 
and XMM-Newton measurements is related to variability 
in the column density, we re-analyzed the Chandra spec- 
trum. In our previous Chandra analysis ( Tomsic k et al.l 



20061). we used the fAnders & Gre vessel ( Il989t) rather than 
Wilms. Allen & McCravl (120001) abundances, and we fitted 



the spectrum with on ly a single power-law. First, w e find 
that when we use the IWilms. AUen & McCravl ( l2000l) abun- 
dances and re-fit the Chandra spectrum, we obtain a value of 
A^H = (5.4;I;j ■') X 10^^ cm"^, which is already somewhat higher 
than found in the previous analysis. Second, the quality of 
the Chandra spectrum does not allow for the detection of 
the soft excess, but we re-fitted the Chandra spectrum after 
including a soft excess with the values of A^h and F mea- 
sured with XMM-Newton (see Table [T]l- This fit gives values 
for the primary power-law component of F = 0.7 ± 0.7 and 
A^H = 2 ^ 10^^ cm"^, which are both consistent with the 
values measured with XMM-Newton. 

Thus, the Chandra spectrum allows for the possibility that 
IGR J 16207-5 129 is an obscured HMXB while the XMM- 
Newton spectrum constrains A^h to be significantly higher 
than the Galactic value, requiring that the source is an ob- 
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scured HMXB. With XMM-Newton, we see that the A^h can 
change, with the first 3,000 s of the observation having A^h = 
(6.4 ± 0.4) X 10^2 cm-2, while the A^h was about a factor 
of two higher for the rest of the observation. The Chan- 
dra spectrum is consistent with either level. Even at the 
highest A^H seen by XMM-Newton, the amount of local ab- 
sorbing material is certainly not as much as seen in some 
of the most extreme systems like IGR J16318-4848, which 
has A^H = 2 X 10^'* cm"^. This is consistent with optical and 
IR observations which have shown P Cygni profiles, for- 
bidden emission lines (indicating a supergiant B[e] spectral 



type) 
4848 


, and IR excesses from local material for IGR J 163 18- 


(iFilliatre & Chatvll2004HMoon et al."2007^.'Rahoui et al. 


2008 


) but not for IGR J16207-5129 dNegueruela & Schurch 


2007 


; iNespoli, Fabresat & Mennickenri2008l; iRahoui et al. 


2008 


). 



4.2. Comparison to other HMXBs Lacking Pulsations 

With the results of our timing study, IGR J 16207-5 129 
joins a relatively short list of accreting HMXBs that show very 
hard energy spectra, resembling known neutron star HMXBs, 
but that do not seem to exhibit pulsations (at least in the 
expected frequency range). One example is 4U 1700-377, 
which has an extremely massive and luminous 06.5 laf-n com- 
panion in a 3.41 day orbit with the compact object. The 
properties of this source as seen by XMM- Newton are remark- 
ably s imilar to those of IGR J 1 6207-5 1 29 (Ivan der Meer et all 
l2005h . Light curves with 10 s time resolution show flares and 
dips where the flux changes by factors of at least 5 on ks time 
scales. While the spectral analysis was limited to the fainter 
time periods because of photon pile-up, spectral fits show a di- 
rect component with a hard power-law index (F = 1.08-1.87 
for the "low-flux" spectrum, depending on the exact model 
used) and a high level of absorption (A^h = 6.8 x 10^^ cm"^ to 
A^H = 2.0 X 10 cm"^) along with a soft excess. While the light 
curves for 4U 1700-377 and IGR J 16207-5 129 are somewhat 
unusual in their extreme variability on long time scales, it 
should be noted that many pulsating HMXBs (e.g.. Vela X-1 
and GX 301-2) also show very similar energy spectra to those 
described here, with hard primary power-law components and 
soft excesses (Nagase 1989). 

A comparison of the power spectra of 4U 1700-377 and 
IGR J 16207-5 129 show similarities also. The typical power 
spectrum of an HMXB pulsar has 3 parts: A relatively flat 
portion at frequencies below the pulsation frequency; a re- 
gion where the pulsations, harmonics, and sometimes quasi- 
periodic oscillations (QPOs) are observed, and a higher fre- 
quency region where the power spectrum can be described 
as a power-law with a slope of 1.4-2.0 (Belloni & Hasingeg 
Il990l: iHoshino & TakeshimallT99l "ciiakrabarty e t al.ii200llK 
However, a uniform study of 12 HMXBs with £X0&4r indi- 
cated that4U 1700-377 (and, interestingly, Cyg X-3) deviated 
from this pattern by showing only a steep power-law in the 
0.002-1 Hz frequen cy range ([B elloni & Hasinger 1990), like 
IGR J16207-5129. Belloni & Hasinger. (,1990) quote 0.008- 
29 Hz integrated rms noise levels for 4U 1700-377 between 
6 and 12%, and if we recalculate the IGR J16207-5129 noise 
level for this frequency range, we obtain 1 1.7%±2.7%. 

While there are many similarities between these two 
sources, one possible difference is that the average lu- 
minosity of 4U 1700-377 is probably higher. To avoid 
model-dependent flux measurements, we compare the aver- 
age 20-40 ke V flux measureme nts made by INTEGRAL. For 
4U 1 700-377, iMetalJ (EQOTI) quote a value of 208.1 ±0.1 



millicrab compared to 3.3 ±0.1 millicrab for IGR J 16207- 
5129. The distanc e to 4U 1700-377 is estimated at 1.9 kpc 
dAnkav et al.ll200 lV so the ratio of fluxes for the two sources 
would imply a distance of ^^15 kpc for IGR J 16207-5 129 
if their 20^0 keV luminosities were the same, which is 
at the very upper limit of the distance range derived by 
INespoli. Fabregat & MennickentI (l2008b . At the best esti- 
mate for die IGR J 16207-5 129 distance, 6.1 kpc, the lumi- 
nosity of IGR J 16207-5 129 would be about 6 times lower 
than 4U 1700-377. At a distance of 6.1 kpc, the aver- 
age IGR J 16207-5 129 20-40 keV luminosity is 1.1 x lO''^ 
ergs s~', and the average 0.5-10 keV unabsorbed luminosity 
is 1.6 x lO-'^ ergs s~'. 

Despite soft X-ray spectra that are similar to neutron star 
HMXBs, the nature of the compact object is still unclear 
for both 4U 1700-377 and IGR J 16207-5 129 because of 
the lack of X-ray pulsations. One reason that the black 
hole possibility has been taken seriously for 4U 1700-377 
is that the compact object mass ha s been measured to be 
2.44 ±0.27 Mq (lClark etalJl2002h . which is significantly 
higher tha n the values close to 1 .4 Mfri measured for most neu- 
tron stars dThorsett & Chakrab artv 1999), and could indicate 
that 4U 1700-377 harbors a low mass black hole. Obtain- 
ing a compact object mass measurement for IGR J16207- 
5129 would certainly be interesting and may be feasible as 
it is relatively bright in the optical and IR (see §1). However, 
the mass measurement would be challenging because accurate 
spectroscopy would be required to measure the massive B 1 la 
companion's radial velocity curve. Rather than black holes, 
these sources may harbor very slowly rotating neutron stars. 
Several such X-ray binarie s are known: IG R J16358-4726 
with its 1.64 hr pulsations dPatel et al.ll2004 : 28 0114+650 
with its 2.7 hr pulsations and B 1 la spectral type (i Farrell et alJ 
l2008l and references therein), which is the same spectral type 
as IGR J16207-5129; 4U 19 54+319 with its Hkely 5 hr pul- 
sations (iMattana et al.l l2006l) ^nd IE 161348-5055.1 with 
its likely 6.7 hr pulsations dPe Luca et al.l l2006). With our 
-12 hr XMM-Newton observation of IGR J16207-5129, we 
can rule out pulsation periods in the —1-2 hr range but prob- 
ably not in the <; 4 hr range. 

Another HMXB with a neutron star-like spectrum that 
has been well-studied without detecting pulsations is 
4U 2206+5 4. This syste m has an 09.5 V or 09.5 III com- 
panion (Rib o etalJl2006l) and a 9.6 or 1 9.25 day orbital pe- 
riod (iCorbet. Markwardt & TueUeiil2007h . Like IGR J 1 6207- 
5129, its X-ray emission is highly variable dMasetti et alj 
l2004l) . but its energy spectrum is significantly different, show- 
ing no evidence for local absorption. Its 10~^-1 Hz power 
spectrum is dominated by s trong red (i.e., power - law) noise 
dNegueruela & Reigl 12001 iTorrejon et aP 120041: iBlay et alJ 
HOOI, which is similar to IGR J 16207-5 129 and 4U 1700- 
377. However, for 4U 2206+54, the possible detection of a 
cyclotron line at ~30 keV provides som e additional evidence 
that the compact object is a neutron star dTorrejon et al.ll2004t 
iBlay et al]l2005l) . If confirmed, this would indicate that the 
lack of pulsations and the red noise power spectrum should 
not be taken as evidence for a black hole. Rather, the non- 
detections of pulsations could be related to a system geome- 
try where, e.g., the neutron star magnetic and spin axes are 
nearly aligned, or the non-detections could be due to spin pe- 
riods that are longer than the duration of the observations. 

Finally, it is worth noting that several IGR HMXBs be- 
sides IGR J 16207-5 129 have not yet shown X-ray pulsations. 
IGR J 1 9 140+095 1 has similarities to IGR J 1 6207-5 1 29, with 
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a B0.5 I spectral type, a high level of obscuration (although 
the level depen ds on orbital phase), and a soft excess (see 
iPrat etalj|2008l and references therein). IGR J16318^848 
(see §4.1) also has not exhibited pulsations despite a high 
level of X-ray variability that has been seen in several soft 
X-ray observations. Although the apparent lack of pulsations 
may be due to the fact that the timing properties of these 
sources have not been very well-studied yet, it is possible that 
the IGR HMXBs include a population of systems with black 
holes or very slowly rotating neutron stars. 

4.3. Soft Excess and Iron Line 

The presence of a soft excess and emission lines are 
very common in both pulsating and non-pulsating obscured 
HMXBs. In addition to iron lines, some sources (e.g.. 
Vela X-1 and 4U 1700-377) have emission lines from Si, 
Ne, etc. in the .5-3 keV range mald i ig up at least part of 
the soft excess (Watanab e et alJ l2006t iBoroson et"a!Tl2003l: 
l^n der Meer et al. 2005). Although our spectral analysis 
does not allow us to definitively rule out other models be- 
sides a power-law (e.g., Bremsstrahlung or black-body) for 
the soft excess, the fact using a power-law leads to a power- 
law index that is consistent with the value of F for the pri- 
mary power-law components allows for two possible physi- 
cal interpretations. One possibility is the partial covering of 
the power-law source by absorbing material (likely the stellar 
wind), which is the scenario envisioned for the spectral fits 
described in §3.2.2. Another possible explanation is that the 
soft excess emission originates as part of the primary power- 
law component, but it is scattered in the stellar wind. In this 
picture, the fact that the soft excess has a much lower column 
density than the primary component would indicate that the 
photons that are part of the soft excess come from the edge 
(i.e., within one optical depth of the edge) of the stellar wind. 

Although it is not clear whether any of the IGR J 16207- 
5129 soft excess comes from emission lines, the fact that the 
flux of the Fe Ka line at 6.4 keV correlates with the overall 
source flux (i.e., the equivalent width is consistent with being 
constant) is different from Vela X-1 and 4U 1700-377, which 
both show lines with very large equivalent widths at low lumi- 
nosities. Regardless of the presence (or not) of lower energy 
lines in IGR J 16207-5 129, the energy of the iron line indi- 
cates a low ionization state implies that the line originates in 
cool material. One possible location for the cool material is 
an accretion disk around the IGR J 16207-5 129 compact ob- 
ject; however, many HMXBs with strong stellar winds have 
compact objects that accrete directly from the wind with lit- 
tle or no disk accretion. The equivalent width of 42 ± 12 eV 
that we measure for the IGR J 16207-5 129 iron line is con- 
sistent with simulations in which a spherical distribution of 
cold matter with solar abundances and A^h ~ 10^^ cm"^ is il- 
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luminated by a central X-ray source ( lMattl | 2 002 [). and this i s 
also ti-ue for several other IGR HMXBs (W alter et alJ l2006). 
Although it is not obvious that the matter in a stellar wind 
will be cold (i.e., neutral), it has been suggested that the wind 
might be a clumpy two-phase medium with cool, dense re- 
gions that are responsible for the emission lines alo ng with 
hot, highly-ionized regions dvan der Meer et al.l2005L and ref- 
erences therein). The partial covering model that we use to 
fit the IGR J 16207-5 129 continuum would also be consistent 
with absorption of the X-ray source by a clumpy wind. 

5. SUMMARY AND CONCLUSIONS 

The XMM-Newton observation of IGR J 16207-5 129 con- 
firms that it is a member of the group of obscured IGR 
HMXBs. We measure a column density of A^h = (1 • 19!j}o5) ^ 
10^^ cm"^ for the average spectrum. We find that the column 
density could have been that high during the only previous 
soft X-ray observation of this source with Chandra, but a de- 
tailed spectral analysis of the XMM-Newton data shows that 
A^H can vary by a factor of ~2. We detect an iron line in 
the energy spectrum, and its strength is consistent with what 
is expected for a spherical distribution of material with the 
measured A^h around the compact object. Although the X- 
ray spectrum is similar to those seen from other neutron star 
HMXBs with a hard primary power-law component and a soft 
excess, we do not detect the pulsations that might be expected 
if the compact object is a neutron star 

A detailed comparison between IGR J 16207-5 129 and an- 
other apparently non-pulsating HMXB, 4U 1700-377, shows 
strong similarities in spectral and timing properties. Most no- 
tably, the power spectra of both sources can be described as 
a single power-law down to 10~^^-10~^ Hz. Since pulsating 
HMXBs show power spectra that break near the pulsation fre- 
quency, it is possible that both of these sources harbor very 
slowly rotating neutron stars (although the possibility that the 
compact object is a black hole cannot be ruled out entirely). 
We note that several of the IGR HMXBs are either known 
to harbor slowly rotating neutron stars or may harbor slowly 
rotating neutron stars in cases where pulsations have not been 
detected. It is possible that in addition to uncovering obscured 
HMXBs, the HMXBs that INTEGRAL is uncovering tend to 
contain slow rotators. 
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TABLE 1 
Spectral Results 



Parameter Value" 

Primary Power-Law Component 

Nn (1.19::«;0«)x 1023cm-2 



FpL (unabsorbed, 0.5-10 keV) (3.68±0.10) x 10''' ergs cm"-^ s'' 



Soft Excess 




1.7 X 10^2 cm-- (fixed) 


r 




FpL 




Iron 


Line Parameters 




1.7 X 10^2 cm-2 (fixed) 


^line 


6.39 ± 0.03 keV 


Cline 


<0.12keV 


N]me 


X Ur^ photons cm'^ s"' 


Equivalent Width* 


42±12eV 




1402/1300 



"The errors on the parameters are for Ax^ = 2.7, corresponding to 90% confidence for one parameters of interest. 

''This is the equivalent width for the unabsorbed Gaussian Une relative to the unabsorbed power-law continuum components. 



TABLE 2 

Spectral Results with the Partial Covering Model vs. Intensity 





Value 


Value 


Value 


Value 


Value 


Parameter 


(aU rates) 


(4-8 c/s) 


(2-4 c/s) 


(1-2 c/s) 


(0-1 c/s) 


Partially Absorbed Power-law Parameters 


Nil (xl023 cm-2) 


1.00^:°^ 


0.64 ±0.04 


1.14 ±0.06 


1.31 ±0.09 


1.24±0.13 


r 


1.13 ±0.02 


1.06 ±0.05 


1.15±0.06 


1.21 ±0.08 


1.30±0.12 


f" 

Fpl'' 
FppL - Fptf^ 
FsE=FpL(l-f)'' 


0.963 ±0.003 

■'•''^-0.04 
-^•'"-0.04 

0.142 ±0.012 


0.964 ±0.008 
11.2±0.3 
10.8 ±0.3 
0.40 ±0.09 


0.980 ± 0.004 0.971 ±0.005 

eoe+o-^'* 3 33+0'* 

"•""-0.11 -'•^-'-0.17 

5.94:::0-2| 3.23 ±0.17 
0.12± 0.02 0.10 ±0.02 


"■^-'^-0.011 
1 9',+0.12 
'•^^-0.10 

1 16+0" 
^•^°-o.io 

0.060 ±0.013 


Iron Line Parameters 


£ii„e (keV) 
(Tiine (keV) 

EW (eV) 


6.39 ±0.03 

06+''"3 

1 6+0-5 
42+13 


6.44 ±0.08 

0.06^ 

4 Q+2.1 

44+18 
^^-30 


6.33 ±0.05 
0.06" 
2.5^0-8 

43+14 


°-^^-0.06 

0.06" 
1.0 ±0.6 

32± 19 


6 39+0- 10 

-0.09 

0.06" 

44!33 




1402/1301 


489/468 


528/503 


432/397 


232/193 



"Fraction of the power-law flux absorbed by interstellar material and material local to the source. A fraction equal to 1-/ is absorbed by just interstellar material 
(at /Vh = 1.7 x 10^^ cm~2), leading to the soft excess seen in the spectrum. 
^Unabsorbed 0.5-10 keV flux in units of 10~^' ergs cm""^ s~^. 

"^Fppi is the unabsorbed 0.5-10 keV flux of the primary power-law component in units of 10~" ergs cm~^ s"' and is derived from Fpi^ and /. 
''Fse is the imabsorbed 0.5-10 keV flux of the soft excess component in units of 10~" ergs cm~^ s~' and is derived from Fpl and /. 
"Fixed. 

Aron Une flux in units of 10"' photons cm~^ s~' . 



